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Thermodynamic quantities are derived for superconducting and pseudogap regimes by taking 
into account both amplitude and phase fluctuations of the pairing field. In the normal (pseudogap) 
state of the underdoped cuprates, two domains have to be distinguished: near the superconducting 
region, phase correlations are important up to the temperature T^. Above T^, the pseudogap region 
is only determined by amplitudes, and phases are uncorrelated. Our calculations show excellent 
quantitative agreement with specific heat and magnetic susceptibility experiments on cuprates. 
We find that the mean field temperature To has a similar doping dependence as the pseudogap 
temperature T*, whereas the pseudogap energy scale is given by the average amplitude above Tc- 



One of the most intriguing problems in high tempera- 
ture superconductivity is the presence of a region above 
the critical temperature Tc and below a temperature T* 
where observable quantities deviate from Fermi hquid be- 
haviour. This region is called pseudogap region be- 
cause it contains effects similar to superconductivity hke 
a partial suppression of electronic density of states. 

The origin of such a pseudogap above Tc is unclear. 
There are four major approaches concerning its theoret- 
ical understanding: the first is based on the formation 
of incoherent Cooper pairs above Tc- Phase order Q or 
Bose condensation would then establish superconduc- 
tivity at Tc- The second assumes that the pseudogap 
is induced by anti-ferromagnetic fluctuations Q. The 
third approach is based on spin-charge separation where 
spins bind together to form spin-singlets and the energy 
needed to split them apart leads to the formation of a 
"spin-gap" [g. The fourth assumes the existence of a 
quantum critical point Q but the latter has never been 
observed. However, these approaches seem to be unable 
to describe specific heat and magnetic susceptibility. 

The main aim of this article is to show that various 
experimental observations can indeed be interpreted in 
terms of fluctuations of the pairing field ijj = and 
that two temperature regions have to be distinguished 
(see Fig. E)): for a relatively small temperature interval 
Tc < T < Tij, the phase of ip is still correlated in space 
over some correlation length ^ (the Kosterlitz-Thouless 
correlation length in 2D). Thus, in this regime, observ- 
ables are governed by correlated phase fluctuations de- 
scribed by the XY- model. For T^ < T < T* , phases of ^ 
are essentially uncorrelated (^ is on the order of the lat- 
tice constant), but \ip\ is still fluctuating and non-zero, 
signaling local pair fluctuations. This explains the wide 
hump seen in specific heat experiments the depres- 
sion of the spin susceptibility Q and the persistence of 
the pseudogap for T < T* . 

Our method has a major difference with respect to the 
Emery and Kivelson phase fluctuations scenario Q of 
the pseudogap regime: we show that phase fluctuations 
influence the pseudogap only up to a temperature T^ 
which is much smaller than T*. Above T^, observables 
are thus only determined by the amplitude of the pairing 
field. 



The picture of two different regimes above Tc is also 
supported by other experiments: a) Demsar et al 0|, 
interpreting the real-time measurements of the quasipar- 
ticle relaxation dynamics, find a temperature interval of 
only a few K in which pair fluctuations associated with 
their collective phase are important, whereas the pseudo- 
gap persists to much higher temperatures, b) Hall effect 
measurements in underdoped GdBa2Cu307_a; show 
a characteristic temperature T' between Tc and T*, at 
which the temperature dependence of cot(OH) deviates 
from T^ and the Hall coefficient has a peak. A possible 
explanation consists in considering vortex excitations as 
scattering centers modifying the Hall angle Qh and Hall 
coefficient, which would again suggest identifying T' with 
our T^ below which correlated vlatex ortices exist. 

We base our calculations on a d-wave attractive Hub- 
bard model 

^ = -*E c^W-t^E^iw QrfW (1) 

{i.j)a i 

with a hopping t between nearest neighbour sites i and j 
on a square lattice. The interaction favours the formation 
of onsite c?-wave pairs since Q\{i) ~ ^ij Olj where 
Dij = 1, (—1) for i being the nearest neighbour site of j in 
horizontal (vertical) direction. Q\- = (c-^c^j^ — c|j^cL)/-\/2 
is the singlet pair operator of neighbouring sites. Decou- 
pling the interaction with the help of a Stratonovich- 
Hubbard transformation, the partition function Z = 
Tr is then 



Z = D^^j(Te-^o'i-^.{m'+^ Qi«+hc)^^ 

where ij) = ipii^r) = |'0(i, T)|e"'^^*''^\ and 7J„ is the non- 
interacting part. The trace over the fermionic operators 
can be evaluated yielding 



Z : 



(2) 



Here G is a Nambu matrix of one-electron Green func- 
tions for fermions interacting with a given, space and 
time dependent pairing field '0(i, r). 

Expanding (|2Jl in powers of V-i/;, Z can be written as 
a functional integral involving an action S[%Ij\ for a field 
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ip that changes slowly in space and that can be taken 
time- independent : 



SM) + 5i(VV') 



(3) 



where 5*0 is a local function of |^'(r)|, and 5*1 = c|V'!/'P/2 
where c is a constant, d is the dimension. 

Now we would like to compute observables such as 
energy, specific heat and spin susceptibility. Our main 
strategy is to treat separately amplitude and phase fluc- 
tuations. In this spirit, two different approaches are pos- 
sible: 1) The amplitude is fixed and determined by a 
suitable variational equation. 2) The energy is expanded 
around the average amplitude. 

Variational Method: we can neglect amplitude cor- 
relations since simulations show that they are weak be- 
tween different sites (iV'illV'jl) — (IV'P) ~ (i-C- the 
amplitude is always positive and shows no critical be- 
haviour) . Rewriting the free energy F in terms of a con- 
stant amplitude yields 
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where the Jacobian of the polar transformation is put 
into the exponential, and V is the volume. Equating the 
derivative of F with respect to |-0| to zero leads to the 
self-consistent equation: 
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The first term, called the amplitude contribution, leads 
to the BCS gap equation if other contributions are 
neglected. The second comes from the Jacobian and im- 
plies that the amplitude is never zero. The third term is 
the expectation value of the energy in the XY model de- 
pending on the constant coupling c\ip\'^ /2. This contribu- 
tion characterizes the influence of the phase fluctuations. 
Solutions of equation ^ are reliable for all temperatures 
except for T << Tc- They are accurate at Tc only in the 
underdoped regime. 

Average Value Method: now, we keep the local cou- 
pling in 5*1 between amplitude and phase by introducing 
the statistical Ginzburg-Landau (GL) model (see [T^). 
Averages like (I?/"]) and (5*1) are computed within the 
GL model with all fluctuations and correlations. Then 
fermionic observables are expanded with respect to the 
average amplitude. Expanding the energy density ob- 
tained from action Q around the average amplitude (see 
Fig. yields: 

E = {S)s/V « Soim) + {Si)s + OiiSm') (5) 

where the square of amplitude fluctuations is neglected 
in first approximation. Averages are computed using a 
normalised GL action: 



SglIiP] = keVo j d\ (Ugl + |V^|V2) 



(6) 
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FIG. 1: Average amplitude (thick) and the correspond- 
ing standard deviation (dV'P) — (long dashed) from 
2D GL simulations. Parameters are Vo = 95K, To = 260K. 
Dashed curves: amplitude distribution P(|7/'|) (histogram) for 
T = 300K (right) and T = 30K (left). 

where Uql = V^[{T/To - l)|i/i(r)|^ + \ip{r)f /2] comes 
from the expansion of Sq in powers /3|'!/'|- To is the mean 
field pairing temperature, and ip = ip/ \ipo\ where ipo = 
Tp{T = 0). S'gl is normalised with a lattice spacing e. 
T] :~ e/^o, where is the mean field correlation length 
at zero temperature. Vq is the zero temperature phase 
stiffness. The energy becomes 



E^Soim)GL) + {Si) 



GL 



(7) 



where 5o((|^|)) = ((|^|>VC^) - ^ log[2cosh/3i?,/2] 
corresponds the BGS free energy for which the gap value 
is determined by the GL average. The d-wave quasi- 
particle energy is Eg = [{sq — /i)^ + cos^(26')]^/^ 
where fj, is the chemical potential, and 6 is the angle in 
k space with respect to direction. 

For computer simulations of the statistical ensemble 
{ijj} under the action Sgl, we use a standard Monte 
Carlo procedure to update amplitude \ip\ and a Wolff 
algorithm for the phase as for the real model [lj|. 
Typically 10* sweeps are needed to obtain good statistics 
for lattice size: N = 40^. The value of rj depends on the 
corse-graining procedure and is fixed to 3 in this method. 
Changes in rj modify slightly the shape of the average am- 
plitude. The link between Sq and the GL action is made 
by fixing ratio of jV'ol Stiid Tq. Since we are interested in 
the temperature domain from Tc to T*, the ratio |'i/'o|/2o 
is fixed at T ~ Tc by requiring that the specific heat 
jump recovers its BCS value in the overdoped regime. 
We found that this ratio must be l^/jol/Tb = 3.06 in the 
s-wave case, and iV'ol/To = 3.72 for c?-wave symmetry. 

Both approaches are valid below and above the critical 
temperature Tc which is the temperature where the phase 
stiffness becomes zero. 

The specific heat C is the sum of Cq and Ci, resp. 
amplitude and gradient contributions. Defining 7 = 
C / {"fnT) where 7„ the Sommerfeld constant, the reduced 
specific heat is 



7 = 7o((|'0l))+7i 



(8) 
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where 71 is divided by Tc instead of T since 5*1 is classical 
and does not satisfy the third law of thermodynamics. 
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FIG. 2: The specific heat from the variational method 
(thick) , which is the sum of gradient (dashed) and amplitude 
(dotted-dashed) contributions, reproduces measurements of 
YBa2Cu306.73 (points). Inset: T dependent amplitude 
from equation Q • 
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FIG. 3: Experimental specific heat (points) for YBa2Cu306+a; 
compared to the average value method. 

The amplitude contribution 70 is 1 at high temperature. 
The gradient contribution is normalised as: 



In 



kB C\ 

NkB 



(9) 



where V = N^q, and cl'^'^ /{NkB) is the specific heat per 
number of lattice sites coming from the simulations. Ex- 
periments give 7„ w 26 mJ K"-'^ mol~^ = 252 J K^-'^ m""^. 
For the fit of Fig. |2 using the reasonable value ^0 ~ 16 A, 
we get the dimensionless constant a = kB / (Sl^^n) ~ 13.5. 
In Fig. |21the experimental specific heat of YBa2Cu306.73 
is fitted using the variational method reproducing 
the double peak structure: a sharp peak at coming 
mainly from correlated phase fluctuations and a wide 
hump rounded by amplitude fluctuations (see Fig. ^ 
for typical distributions). The crossover temperature Ti 



corresponds to the temperature where 71 is less than ap- 
proximately 2% of the normal specific heat. In the am- 
plitude equation Q), a 2-dimensional density of states 
D{e) = l/W is used with W = 5000K, = 0.25VF 
and U = 959K. These values give Tq « 200K and 
■00 ~ 2.1ATf) in agreement with experiments The 
other parameters are Vq = 72K and rj = 5. 

The average value method is used to reproduce specific 
heat measured for different dopings in Fig. 13 For un- 
derdoped systems x < 0.80 we use simulations in D = 2. 
For the more overdoped, x > 0.80, simulations are done 
in D = 3. Parameters Vq and Tq extracted from the 
best fits are shown in the phase diagram of Fig. [3 a 
ranges from 7 to 15. The magnetic susceptibility x 
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FIG. 4: The measured spin susceptibility of YBa2Cu306.63 
(points) divided by xo is well fitted by the theoretical (thick) 
susceptibility Xv- The dotted-dashed line is standard devia- 
tion of the average amplitude. 

has two contributions: the paramagnetic spin suscepti- 
bility Xv orbital diamagnetic susceptibility Xd- 
X-p has been measured by Takigawa et al ^ on powder 
YBa2Cu306.63 using Cu and O NMR experiments. The 
direct contribution coming from phases is negligible be- 
cause NMR probes essentially the presence of pairs which 
is related to the amplitude. Therefore. Xp is given by the 
amplitude contribution: 



Xp 



2T 



2^ 



de cosh 



-2 



(|^|)g.^cos2(20) 



2T 



(10) 

Xo is the Pauli spin susceptibility. In Fig. 21 we com- 
pare the result of equation l|l()() and Takigawa's mea- 
surements on powder YBa2Cu3 06.63- The best fit yields 
To = 159. 2K, Vb = 59.3K. The competition between 
amplitude and thermal energy enters into the spin sus- 
ceptibility Hl()|l by the ratio {\il}\)/T. The temperature 
where phases start to influence the thermodynamics 
is the temperature where deviates from the average 
amplitude computed for random phases. This tempera- 
ture is found above Tc at T^ ~ 90K. The orbital dia- 
magnetic susceptibility Xd in YBa2Cu3O6.60 shows fluc- 
tuations effects up to about 15 K above T^- For T > T^, 
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phase fluctuations arc so strong that Xd vanishes jl^ . 
The Green function can be calculated in CPA (Co- 
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FIG. 5: Phase diagram of YBa2Cu306+a:. Effects of am- 
plitudes are large in the quasi uncorrelated pseudogap region 
(below To), whereas phase correlations remain important only 
below r^. The temperature T* , where 7 and Xd cross over to 
normal behaviour, is located in the hatched area. Inset: the 
pseudogap energy scale Eg. 

herent Potential Approximation) for electrons that are 
scattered from spatially uncorrelated " pair impurities" . 
We extended the s-wave CPA [13 by using a d-wave am- 
plitude and uncorrelated phases. We have computed the 
density of states at half-filling fi = 0.514^ and maximum 
amplitude |?/'| = 0.81^. The result is a pseudogap of 
width A « 0.2iy although phases are completely uncor- 
related. 

Phase diagram: Values of To, Vo and extracted from 
specific heat fits (Fig. O arc reported in Fig. [S] Ti is the 
temperature where Ci/7„ = 2%. is computed in the 
same way as shown in Fig. 0] using average amplitudes 



from specific heat fits. Here, it is the temperature where 
the amplitude differs by 2% from the amplitude in a ran- 
dom phase field. The energy scale Eg := {\iP\)t=200K of 
the pseudogap is defined as the amplitude at T = 200K. 
Eg shows the same doping dependence as the one found 
by Loram and Tallon |l9j . However it is not related to 
any hidden critical point. Phase correlations above Tc 
grow rapidcly in the underdoped regime following the 
line, an have a similar doping dependence as Nernst 
effect results 0| . The gradient specific heat from Si dis- 
appears more rapidcly like in the Hall effect ^3 ■ 
Discussion: amplitude and phase fluctuations are the 
key for understanding the pseudogap regime of under- 
doped high temperature superconductors. Phase correla- 
tions disappear completely near a temperature above 
Tc, and therefore, for T > T,f,, the pseudogap region is 
dominated by amplitude fluctuations. The mean field 
temperature Tq has a similar doping dependence as T*, 
signaling that the pseudogap region is due to indepen- 
dent fluctuating pairs. Comparisons with measured spe- 
cific heat on underdoped YBCO reproduce the double 
peak structure: a sharp peak at Tc coming mainly from 
phase correlations and a separate wide hump below T* 
rounded by fluctuations. The spin susceptibility, related 
to the amplitude, recovers its normal behaviour near T* 
whereas the orbital magnetic susceptibility, related to 
phases, disappears near T^. These considerations are in- 
dependent from the underlying pairing mechanism, and 
any microscopic theory inducing pairing should lead to 
similar conclusions. 

All these findings provide additional evidence for the fact 
that superconductivity and pseudogap have the same ori- 
gin. The former is primarily related to phases of the 
pairing field, which are ordered below the transition tem- 
perature and whose correlations survive over a limited 
temperature region above Tc- The pseudogap regime of 
underdoped materials then extends to much higher tem- 
peratures due to the persisting amplitude fluctuations of 
the pairing field. 
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